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1 Introduction 

A classic question in probability theory, studied by M. Kac, S. O. Rice and 
many others, is to find the expected number and distribution of zeroes or 
critical points of a random polynomial. The same questions can be asked for 
random holomorphic functions or sections of bundles, and are the subject of 
"random algebraic geometry." 

While this theory has many physical applications, in this article we focus on a 
variation on a standard question in the theory of disordered systems. This is 
to find the expected distribution of minima of a potential function randomly 
chosen from an ensemble, which might be chosen to model a crystal with 
impurities, a spin glass, or another disordered system. Now whereas standard 
potentials are real-valued functions, analogous functions in supersymmetric 
theories, such as the superpotential and the central charge, are holomorphic 
sections of a line bundle. Thus one is interested in finding the distribution 
of critical points of a randomly chosen holomorphic section. 

Two related and much-studied problems of this type are the problem of 
finding attractor points in the sense of Ferrara, Kallosh and Strominger, 
and the problem of finding flux vacua as posed by Giddings, Kachru and 
Polchinski. These problems involve a good deal of fascinating mathematics 
and are good illustrations of the general theory. 
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2 Elementary random algebraic geometry 



Let us introduce this subject with the problem of finding the expected dis- 
tribution of zeroes of a random polynomial, 



f(z) = C + CiZ + . . . c N z N . 

We define a random polynomial to be a probability measure on a space of 
polynomials. A natural choice might be independent Gaussian measures on 
the coefficients, 

N 

dM = dv[c , , c N ] = n d 2 Cl |I e -W W (1) 

i=0 

We still need to choose the variances. At first the most natural choice would 
seem to be equal variance for each coefficient, say = 1/2. We can charac- 
terize this ensemble by its two-point function, 

G( Zl ,z 2 ) = E[f( Zl )r(z 2 )} = ]>[/] fWfr) 

1-Z1Z2 

We now define d[/,o(z) to be a measure with unit weight at each solution 
of f(z) = 0, such that its integral over a region in C counts the expected 
number of zeroes in that region. It can be written in terms of the standard 
Dirac delta function, by multiplication by a Jacobian factor, 

dto(z)=E[6<V(f(z))df(z)dr(z)]. (2) 

To compute this expectation value, we introduce a constrained two-point 
function, 

E[sw(f(z)) m r( Z2 )} 

G f(z)=0 {z 1: z 2 ) E[6W{f(z))] 

It could be explicitly computed by using the constraint f(z) =0 to solve for 
a coefficient q in the Gaussian integral, i.e. projecting on the linear subspace 
= ^CiZ 1 . The result, in terms of G(z 1: z 2 ), is 
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, _ _ . _ G(zi, z)G(z, z 2 , 

Uf^ z ) =0 [Z 1 , Z 2 ) — Lr{Z 1 , Z 2 ) — ■ 



G(z,z) 

as can be verified by considering 

E[5^(f(z)) f(z) r(z 2 )) oc G m=0 (z,z 2 ) = G(z,z 2 ) - ~ = 

Eq. (J2J) follows from this simply by taking derivatives: 

^oO) = ccb) ^ m z llZ 2^z D 1 D 2 G z (z 1 , z 2 ) 
= ±dd log G{z,z). 

For the constant variance ensemble Eq. (J2J), 

, , . d 2 z ( 1 (N + l) 2 (zz) N \ 
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We see that as A" — > oo, the zeroes concentrate on the unit circle \z\ = 1 
(Hammersley, 1954). 

A similar formula can be derived for the distribution of roots of a real poly- 
nomial on the real axis, using dfi(t) = K[5(f (t))\df / dt\}. One obtains (Kac, 
1943): 
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t 2 ) 2 (l-t 2N+2 ) 2 ' 



Integrating, one finds the expected number of real zeroes of a degree A" 
random real polynomial is ~ -logA^, and as A" — > oo the zeroes are 
concentrated at t — ±1. 

While concentration of measure is a fairly generic property for random poly- 
nomials, it is by no means universal. Let us consider another Gaussian 
ensemble, with variance a n = N\/n\(N — n)\. This choice leads to a partic- 
ularly simple two-point function, 

G(z,z) = (l+zz) N , (4) 

and the distribution of zeroes 



Rather than concentrate the zeroes, in this ensemble zeroes are uniformly 
distributed according to the volume of the Fubini-Study (S77(2)-invariant) 
Kahler metric 

uj = ddK- K = log(l + zz) 
on complex projective space CP 1 . 

We can better understand the different behaviors in our two examples by 
focusing on a hermitian inner product (/, g) on function space, associated to 
the measure Eq. (0) by the formal expression 

dM = [Df]e-v>n. 

In making this precise, let us generalize a bit further and allow / to be a 
holomorphic section of a line bundle £, say O(N) over CP 1 in our examples. 
We then choose an orthonormal basis of sections (sj, Sj) = 5ij, and write 

f = CiSi ( 6 ) 

i 

and 

1 N 

We can then compute the two-point function 

N 

G( Zl ,z 2 ) = E[s( Zl )s*(z 2 )} =Y,Si(z 1 )s* i (z 2 ). (7) 

i=l 

and proceed as before. 

In these terms, the simplest way to describe the measure for our first example 
is that it follows from the inner product on the unit circle, 

Thus we might suspect that this has something to do with the concentration 
of Eq. (J2J) on the unit circle. Indeed, this idea is made precise and generalized 
in (Shiffman and Zelditch, 2003). 

Our second example belongs to a class of problems in which M. is compact 
and £ positive. In this case, the space H°(JA,£) of holomorphic sections is 
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finite dimensional, so we can take the basis to consist of all sections. Then, 
if Ai is in addition Kahler, we can derive all the other data from a choice of 
hermitian metric h(f,g) on C. In particular, this determines a Kahler form 
uj as the curvature of the metric compatible connection, and thus a volume 
form VoLj = uj n /n\. We then define the inner product to be 



Thus, the measure Eq. and the final distribution Eq. (j2J) are entirely 
determined by h. In these terms, the underlying reason for the simplicity 
of Eq. (jSJ) is that we started with the SU(2) invariant metric h, so the final 
distribution must be invariant as well. More generally, Eq. (JJJ) is a Szego 
kernel. Taking C = Cf N for N large, this has a known asymptotic expansion, 
enabling a rather complete treatment (Zelditch, 2001). 

Our two examples also make the larger point that a wide variety of distribu- 
tions are possible. Thus we must put in some information about the ensemble 
of random polynomials or sections which appear in the problem at hand, to 
get convincing results. 

The basic computation we just discussed can be vastly generalized: to mul- 
tiple variables, multipoint correlation functions, many different ensembles, 
and different counting problems. We will discuss the distribution of critical 
points of holomorphic sections below. 

3 The attractor problem 

We now turn to our physical problems. Both are posed in the context of 
compactification of the type lib superstring theory on a Calabi-Yau three- 
fold M. This leads to a four dimensional effective field theory with N = 2 
super symmetry, determined by the geometry of M. 

Let us begin by stating the attractor problem mathematically, and afterwards 
give its physical background. We begin by reviewing a bit of the theory of 
Calabi-Yau manifolds. By Yau's proof of the Calabi conjecture, the moduli 
space of Ricci-flat metrics on M is determined by a choice of complex struc- 
ture on M, denote this J, and a choice of Kahler class. Using deformation 
theory, it can be shown that the moduli space of complex structures, denote 
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this Ai c (M), is locally a complex manifold of dimension h 2,1 (M). A point J 
in Ai c (M) picks out a holomorphic three-form Qj G H 3 '°(M, C), unique up 
to an overall choice of normalization. The converse is also true; this can be 
made precise by defining the period map Ai c (M) — > F(H 3 (M, Z) <g)C) to be 
the class of f2 in H 3 (M, Z) <g> C up to projective equivalence. One can prove 
that the period map is injective (the Torelli theorem), locally in general and 
globally in certain cases such as the quintic in CV 4 . 

Now, the data for the attractor problem is a charge, a class 7 G H 3 (M, Z). 
An attractor point for 7 is then a complex structure J on M such that 



This amounts to h 2,1 complex conditions on the h 2,1 complex structure mod- 
uli, so picks out isolated points in Ai c (M), the attractor points. 

There are many mathematical and physical questions one can ask about at- 
tractor points, and it would be very interesting to have a general method to 
find them. As emphasized by G. Moore, this is one of the simplest problems 
arising from string theory in which integrality (here due to charge quanti- 
zation) plays a central role, and thus it provides a natural point of contact 
between string theory and number theory. For example, one might suspect 
that attractor Calabi-Yau's are arithmetic, i.e. are projective varieties whose 
defining equations live in an algebraic number field. This can be shown to 
always be true for K3 x T 2 , and there are conjectures about when this is 
true more generally (Moore, 2004). 

A simpler problem is to characterize the distribution of attractor points in 
Ai c (M). As these are infinite in number, one must introduce some control 
parameter. While the first idea which might come to mind is to bound the 
magnitude of 7, since the intersection form on if 3 (M, Z) is antisymmetric, 
there is no natural way to do this. A better way to get a finite set is to bound 
the period of 7, and consider the attractor points satisfying 



As an example of the type of result we will discuss below, one can show that 
for large Z max , the density of such attractor points asymptotically approaches 
the Weil-Peterson volume form on Ai c . 



7 G Hj'°{M, C) © Hj'\M, C). 



(8) 



zL x >\z(r,z)\ 



2 _ 



(9) 
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We now briefly review the origins of this problem, in the physics of 1/2 BPS 
black holes in N = 2 supergravity. We begin by introducing local complex 
coordinates z % on M. C (M). Physically, these can be thought of as massless 
complex scalar fields. These sit in vector multiplets of N = 2 supersymmetry, 
so there must be h 2,1 (M) vector potentials to serve as their bosonic partners 
under supersymmetry These appear because the massless modes of the type 
lib string include various higher rank p-form gauge potentials, in particular 
a self-dual four-form which we denote C . Self-duality means that dC = *dC 
up to non-linear terms, where * is the Hodge star operator in ten dimensions. 

Now, Kaluza-Klein reduction of this four-form potential produces b 3 (M) one- 
form vector potentials Aj in four dimensions. Given an explicit basis of three- 
forms uji for H 3 (M, M.) fl H 3 (M, Z), this follows from the decomposition 

h 

C = Aj A uji + massive modes. 
i=i 

However, because of the self-duality relation, only half of these vector po- 
tentials are independent; the other half are determined in terms of them by 
four- dimensional electric-magnetic duality. Explicitly, given the intersection 
form rjij on H 3 ® H 3 , we have 

dAi = rjij * 4 dAj (10) 

where *4 denotes the Hodge star in d — 4. Thus we have h 2,1 + 1 independent 
vector potentials. One of these sits in the N = 2 supergravity multiplet, and 
the rest are the correct number to pair with the complex structure moduli. 

We now consider 1/2 BPS black hole solutions of this four dimensional N = 2 
theory. Choosing any S 2 which surrounds the horizon, we can define the 
charge 7 as the class in H 3 (M, Z) which reproduces the corresponding mag- 
netic charges 

Qi = — dAi = / u>i A 7. 
2tt J S 2 J m 

Using Eq. (JTUJ), this includes all charges. 

One can show that the mass M of any charged object in supergravity satisfies 
a BPS bound, 

M 2 >\z( r ,z)\ 2 . (11) 
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The quantity \Z(j; z)\ 2 , denned in Eq. (P), depends explicitly on 7, and im- 
plicitly on the complex structure moduli z through Q. A 1/2 BPS solution 
by definition saturates this bound. 

We now explain the "attractor paradox." According to Bekenstein and Hawk- 
ing, the entropy of any black hole is proportional to the area of its event 
horizon. This area can be found by finding the black hole as an explicit so- 
lution of four- dimensional supergravity, which clearly depends on the charge 
7. In fact, we must fix boundary conditions for all the fields at infinity, in 
particular the complex structure moduli, to get a particular black hole solu- 
tion. Now, normally varying the boundary conditions varies all the data of 
a solution in a continuous way. On the other hand, if the entropy has any 
microscopic interpretation as the logarithm of the number of quantum states 
of the black hole, one would expect e s to be integrally quantized. Thus, it 
must remain fixed as the boundary conditions on complex structure mod- 
uli are varied, in contradiction with naive expectations for the area of the 
horizon, and seemingly contradicting Bekenstein and Hawking. 

The resolution of this paradox is the attractor mechanism. Let us work 
in coordinates for which the four- dimensional metric takes the form 



With some work, one can see that in the 1/2 BPS case, the equations of 
motion imply that as r decreases, the complex structure moduli z follow 
gradient flow with respect to |Z(7, z)\ 2 in Eq. (fTTj) . and the area A(r) of an 
S 2 at radius r decreases. Finally, at the horizon, z reaches a value at 
which |Z(7, z*)| 2 is a local minimum, and the area of the event horizon is 
A = 47r|Z(7, z*)| 2 . Since z* is determined by minimization, this area will not 
change under small variations of the initial z, resolving the paradox. 

A little algebra shows that the problem of finding non-zero critical points 
of \Z( r y ) z)\ 2 ) is equivalent to that of finding critical points D^Z = of the 
period associated to 7, 



usually called the central charge, with respect to the covariant derivative 



ds 2 = -f{r)dt 2 + dr 2 + 



A(r) 



dttl 




(12) 



D t Z = diZ + (diK)Z. 



(13) 
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Here 

e~ K = JnAti. (14) 

The mathematical significance of this rephrasing is that K is a Kahler po- 
tential for the Weil- Peterson Kahler metric on Ai c (M), with Kahler form 
uj = ddK, and Eq. (|13|) is the unique connection on 

H^°\M,C) regarded 

as a line bundle over Ai c {M), whose curvature is — u. These facts can be 
used to show that DiVt provides a basis for iif( 2,1 )(M, C), so that the critical 
point condition forces the projection of 7 on H^ 2,1 ^ to vanish. This justifies 
our original definition Eq. (jHJ). 



4 Flux vacua in lib string theory 

We will not describe our second problem in as much detail, but just give the 
analogous final formulation. In this problem, a "choice of flux" is a pair of 
elements of H 3 (M, Z), or equivalently a single element 

F G H 3 (M, Z © rZ), (15) 

where r EH = {r G C|Im r > 0} is the so-called "dilaton-axion." 

A flux vacuum is then a choice of complex structure J and r for which 

F G H 3 y°(M, C) © H)' 2 (M, C). (16) 

Now we have h 2,1 + h 0,3 = h 2,1 + 1 complex conditions on the joint choice of 
h 2 ' 1 complex structure moduli and r, so this condition also picks out special 
points, now in M. c x TC. 

The critical point formulation of this problem is that of finding critical points 
of 

W = JttAF (17) 
under the covariant derivatives Eq. (|13|) and 

D T W = d T W + (d T W)Z 

with K the sum of Eq. (|T4^) and the Kahler potential — log Im r for the metric 
on the upper half plane of constant curvature —1. 
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This is a sort of complexified version of the previous problem and arises 
naturally in lib compactification by postulating a non-zero value F for a 
certain three-form gauge field strength, the flux. The quantity Eq. (JT7j) is 
the superpotential of the resulting N — 1 supergravity theory, and it is a 
standard fact in this context that supersymmetric vacua (critical points of 
the effective potential) are critical points of W in the sense we just stated. 

We can again pose the question of finding the distribution of flux vacua in 
M. C {M) x TC. Besides \W\ 2 , which physically is one of the contributions to 
the vacuum energy, we can also use the "length of the flux" 



as a control parameter, and count flux vacua for which L < L max . In fact, 
this parameter arises naturally in the actual lib problem, as the "orientifold 
three-plane charge." 

What makes this problem particularly interesting physically is that it (and its 
analogs in other string theories) may bear on the solution of the cosmological 
constant problem. This begins with Einstein's famous observation that the 
equations of general relativity admit a one parameter generalization, 



Physically, the cosmological constant A is the vacuum energy, which in our 
flux problem takes the form A = . . . — 3|iy| 2 (the other terms are inessential 



Cosmological observations tell us that A is very small, of the same order as 
the energy of matter in the present era, about 10~ 122 Mp lanck in Planck units. 
However, in a generic theory of quantum gravity, including string theory, 
quantum effects are expected to produce a large vacuum energy, a priori of 
order Mp lanck . Finding an explanation for why the theory of our universe is 
in this sense non-generic is the cosmological constant problem. 

One of the standard solutions of this problem is the "anthropic solution," 
initiated in work of Weinberg and others, and discussed in string theory in 
(Bousso and Polchinski, 2000). Suppose that we are discussing a theory with 
a large number of vacuum states, all of which are otherwise candidates to 
describe our universe, but which differ in A. If the number of these vacuum 




(18) 



R»u - ^QiiuR = StcT^ + Ag^. 



for us here). 
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states were sufficiently large, the claim that a few of these states realize a 
small A would not be surprising. But one might still feel a need to explain 
why our universe is a vacuum with small A, and not one of the multitude 
with large A. 

The anthropic argument is that, according to accepted models for early cos- 
mology, if the value of |A| were even 100 times larger than what is observed, 
galaxies and stars could not form. Thus, the known laws of physics guaran- 
tee that we will observe a universe with A within this bound; it is irrelevant 
whether other possible vacuum states "exist" in any sense. 

While such anthropic arguments are controversial, one can avoid them in 
this case by simply asking whether or not any vacuum state fits the observed 
value of A. Given a precise definition of vacuum state, this is a question 
of mathematics. Still, answering it for any given vacuum state is extremely 
difficult, as it would require computing A to 10~ 122 precision. But it is not 
out of reach to argue that out of a large number of vacua, some of them are 
expected to realize small A. For example, if we could show that the number 
of otherwise physically acceptable vacua was larger than 10 122 , and that the 
distribution of A among these was approximately uniform over the range 
(—Mp lanck , Mp lanck ), we would have made a good case for this expectation. 

This style of reasoning can be vastly generalized, and given favorable as- 
sumptions about the number of vacua in a theory, could lead to falsifiable 
predictions independent of any a priori assumptions about the choice of vac- 
uum state (Douglas, 2003). 



5 Asymptotic counting formulae 

We have just defined two classes of physically preferred points in the complex 
structure moduli space of Calabi-Yau threefolds, the attractor points and the 
flux vacua. Both have simple definitions in terms of Hodge structure, Eq. (jHJ) 
and Eq. (JTHJ), and both are also critical points of integral periods of the 
holomorphic three-form. 

This second phrasing of the problem suggests the following language. We 
define a random period of the holomorphic three-form to be the period 
for a randomly chosen cycle in H 3 (M, Z) of the types we just discussed (real 
or complex, and with the appropriate control parameters). We are then 
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interested in the expected distribution of critical points for a random period. 
This brings our problem into the framework of random algebraic geometry. 

Before proceeding to use this framework, let us first point out some differ- 
ences with the toy problems we discussed. First, while Eq. (jTJj) and Eq. (JT7J) 
are sums of the form Eq. (jHJ), we take not an orthonormal basis but instead 
a basis s« of integral periods of Q. Second, the coefficients q are not nor- 
mally distributed but instead drawn from a discrete uniform distribution, i.e. 
correspond to a choice of 7 in H 3 (M, Z) or F as in Eq. (J15)) . satisfying the 
bounds on \Z\ or L. Finally, we do not normalize the distribution (which 
is thus not a probability measure) but instead take each choice with unit 
weight. 

These choices can of course be modified, but are made in order to answer 
the question, how many attractor points (or flux vacua) sit within a specified 
region of moduli space. The answer we will get is a density \±{Z max ) or 
fJ>{L max ) on moduli space, such that as the control parameter becomes large, 
the number of critical points within a region R asymptotes to 



The key observation is that to get such asymptotics, we can start with a 
Gaussian random element of H 3 (M, M) (or flux). In other words, we neglect 
the integral quantization of the charge or flux. Intuitively, this might be 
expected to make little difference in the limit that the charge or flux is large, 
and in fact one can prove that this simplification reproduces the leading large 
L or \Z\ asymptotics for the density of critical points, using standard ideas 
in lattice point counting. 

This justifies starting with a two-point function like Eq. (JJJ). While the 
integral periods Sj of Q can be computed in principle (and have been in many 
examples) by solving a system of linear PDE's, the Picard-Fuchs equations, 
it turns out that one does not need such detailed results. Rather, one can 
use the following ansatz for the two-point function, 




G( Zl ,z 2 ) =Y.i = W J si{z l )s: J {z 2 ) 
= f M n(z 1 )Ati(z 2 ) 
= exp-K(z 1 ,z 2 ). 
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In words, the two-point function is the formal continuation of the Kahler 
potential on Ai c (M) to independent holomorphic and antiholomorphic vari- 
ables. This incorporates the quadratic form appearing in Eq. (j!8|) and can 
be used to count sections with such a bound. 

We can now follow the same strategy as before, by introducing an expected 
density of critical points, 

dfi(z) =E[5 in) {Dis(z))5 in) (Dis(z)) | det %|], (19) 

l<i,j<2n 

where the "complex Hessian" H is the 2n x 2n matrix of second derivatives 

f diDjs{z) diD jS (z)\ . v 

n ~ y^D-fsiz) diD jS {z)) [M) 

(note that dDs = DDs at a critical point). One can then compute this 
density along the same lines. The holomorphy of s implies that diDjS = oufjS, 
which is one simplification. Other geometric simplifications follow from the 
fact that Eq. (|19|) depends only on s and a finite number of its derivatives at 
the point z. 

For the attractor problem, using the identity 

D^s = T ijk uo kk D- k s = 0, 

from special geometry of Calabi-Yau threefolds, the Hessian becomes trivial, 
and det if = \s\ 2n . One thus finds (Denef and Douglas, 2004) that the 
asymptotic density of attractor points with large \Z\ < Z max in a region R is 

on+l 

N{R, \Z\ < Z max ) ~ Z»£ ■ vol(R) 

[n + l)7r n 

where vol(R) = j R u n /n\ is the volume of R in the Weil- Peterson metric. The 
total volume is known to be finite for Calabi-Yau threefold moduli spaces, 
and thus so is the number of attractor points under this bound. 

The flux vacuum problem is complicated by the fact that DDs is non-zero 
and thus the determinant of the Hessian does not take a definite sign, and im- 
plementing the absolute value in Eq. (|T§|) is nontrivial. The result (Douglas, 
Shiffman and Zelditch, 2004) is 

fi(z) 1 / | det(M* - \x\ 2 ■ l)\e HtA ^~ lH -W 2 dHdx 

6 3 ! V det A(z) Jn(z)xc 
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where TC(z) is the subspace of Hessian matrices Eq. (|20p obtainable from 
periods at the point z, and A(z) is a covariance matrix computable from the 
period data. 

A simpler lower bound for the number of solutions can be obtained by instead 
computing the index density 

fjLj(z) =E[6^ n \D i s)S^(D i s) det H tJ }, (21) 

l<i ,j<2n 

so-called because it weighs the vacua with a Morse- Witten sign factor. This 
admits a simple explicit formula (Ashok and Douglas, 2004), 

I vac (R, L < L max ) ~ ^^Tr ■ f det ( n + u ■ *)> ( 22 ) 

where 1Z is the n + 1 x n + 1 dimensional matrix of curvature two-forms for 
the Weil-Peterson metric. 

One might have guessed at this density by the following reasoning. If s 
had been a single- valued section on a compact M. c (it is not), topological 
arguments determine the total index to be [c n+ i(£ <S>T*Ji4)], and this is the 
simplest density constructed solely from the metric and curvatures in the 
same cohomology class. 

It is not in general known whether this integral over Calabi-Yau moduli space 
is finite, though this is true in examples studied so far. One can also control 
\W\ 2 as well as other observables, and one finds that the distribution of \W\ 2 
among flux vacua is to a good approximation uniform. Considering explicit 
examples, the prefactor in Eq. (|2*2*|) is of order 10 100 -10 300 , so assuming that 
this factor dominates the integral, we have justified the Bousso-Polchinski 
solution to the cosmological constant problem in these models. 

The finite L corrections to these formulae can be estimated using van der 
Corput techniques, and are suppressed by better than the naive L -1 / 2 or 
l-Z"! -1 one might have expected. However the asymptotic formulae for the 
numbers of flux vacuum break down in certain limits of moduli space, such 
as the large complex structure limit. This is because Eq. ()18|) is an indefinite 
quadratic form, and the fact that it bounds the number of solutions at all is 
somewhat subtle. These points are discussed at length in (Douglas, Shiffman 
and Zelditch, 2005). 
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Similar results have been obtained for a wide variety of flux vacuum counting 
problems, with constraints on the value of the effective potential at the min- 
imum, on the masses of scalar fields, on scales of supersymmetry breaking, 
and so on. And in principle, this is just the tip of an iceberg, as the study 
of more or less any class of superstring vacua leads to similar questions of 
counting and distribution, less well understood at present. Some of these are 
discussed in (Douglas, 2003; Acharya et al 2005; Denef and Douglas 2005; 
Blumenhagen et al 2005). 
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